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Abstract 

Qualitative features of Ginsparg- Wilson fermions, as formulated by 
Neuberger, coupled to two dimensional U(l) gauge theory are studied. 
The role of the Wilson mass parameter in changing the number of mass- 
less flavors in the theory and its connection with the index of the Dirac 
operator is studied. Although the index of the Dirac operator is not 
related to the geometric definition of the topological charge for strong 
couplings, the two start to agree as soon as one goes to moderately 
weak couplings. This produces the desired singularity in the quenched 
chiral condensate which appears to be very difficult to reproduce with 
staggered fermions. The fermion determinant removes the singularity 
and reproduces the known chiral condensate and the meson mass within 
understandable errors. 
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1. Introduction 



Our understanding of chiral symmetry on the lattice has matured recently. The 
overlap formulation |L] has lead to a new Dirac operator that is expected to describe 
exactly massless fermions 0. The chiral symmetry on the lattice that makes this 
possible has been discovered in Ref. ||. This symmetry arises when the Dirac 
operator D satisfies^ 

D l5 + l5 D = D l5 D, (1) 

called the Ginsparg- Wilson relation j|. This relation appears to be quite powerful. 
It helps to reproduce the continuum relations arising due to chiral symmetry, and 
thus leads to elegant renormalization properties |J . Further, like in the continuum, 
it also suggests an elegant resolution of the U(l) problem on the lattice which has 
interesting consequences for the lattice calculations of QCD ||. We will refer to the 
fermion formulations which obey the Ginsparg- Wilson relation as Ginsparg- Wilson 
fermions. 

In addition to the Neuberger operator [0, the fixed point Dirac operator also 
obeys the Ginsparg- Wilson relation || . Here we do not discuss the latter approach. 
In the context of the Neuberger's Dirac operator it is important to point out that 
there is another closely related approach to massless fermion based on an extension 



of an idea of Kaplan pi called domain wall fermions. These fermions can be 



considered as a truncated version of the overlap (TT| . The appropriate Dirac operator 
of domain wall fermions satisfy the Ginsparg- Wilson relation up to exponentially 
small violations in the limit where the separation of the walls, on which the physical 
fermions live, becomes large. The domain wall fermions have been used in QCD 
calculations Jl2[ and appear to have very good chiral properties compared to the 
Wilson fermion approach. More recently this approach has also been shown to 
reproduce some of the singularities in the chiral condensate due to exact zero modes 
that have not been seen in staggered fermion formulations fl3fl . Thus the domain wall 
fermions appear to be a practical way to implement the Ginsparg- Wilson relation 
approximately. However, in this article we explore the possibility of keeping the 
Ginsparg- Wilson relation exact by working directly with Neuberger's Dirac operator. 

Although chiral symmetry emerges naturally with Ginsparg- Wilson fermions, 
their dynamical consequences are just beginning to be explored. The first study 
of the fixed point Dirac operator was done in the Schwinger model in Ref. |14] 



A similar study for Neuberger's Dirac operator is lacking^. Tests performed in 
closely related formulations like the overlap [16[] and domain wall fermions |L7| can 



only be used to gain confidence, since all these formulations differ in their cutoff 



3 The relation is more general. We only state the simplest form here. 

4 Look at the work by Farchioni et al., pal, which appeared after the completion of our present 



work. 
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effects. In Neuberger's Dirac operator such effects enter through a mass parameter 
M. It is important to understand such effects and check for universality. This has 
not yet been done. Further, in the context of the Schwinger model, the studies 
in the overlap and domain wall formulations have only concentrated on the chiral 
condensate [|18|, [T9|]p| . The meson mass has been ignored. As we discuss below, this 
is another important quantity to calculate to confirm universality. 

A new feature of the Ginsparg- Wilson fermions is that the Dirac operator has a 
non-zero coupling between any two lattice points which are exponentially small in 
the distance between the points |2(J . Such actions are still considered localQ if they 



do not destroy the mass spectrum of the theory. This can be checked in the free 
theory. On the other hand it is important to see if this feature continues at a non- 
perturbative level and no spurious long range correlations originate. This question 



has been recently addressed for smooth gauge fields in Ref. pl| . Non-perturbatively 
one can study the mass spectrum of the theory. Thus determining the meson mass 
in the Schwinger model is important to confirm our beliefs. 

In this article we study a two dimensional U(l) lattice gauge theory coupled to 
fermions proposed by Neuberger ||. In the continuum limit, this model approaches 
the Schwinger model up to four-Fermi couplings 11191^ To control these couplings 
is quite difficult, and will not be attempted here. Instead we will try to learn 
qualitative features of the model as a way to learn how it could reproduce the 
main properties of the Schwinger model if a more thorough and complete analysis 
is made. We also leave the more interesting multi- flavor case to later studies since 
there the dynamics is considerably more complicated and finite volume effects must 
be carefully studied before the effects of complete chiral symmetry guaranteed by 
Neuberger's Dirac operator can be understood. In particular it is important to 
obtain results like (ipip) ~ m 1 / 3 and the pion mass ~ m 2 / 3 [p2|, |17|j . 

This article is organized as follows. In section 2 we introduce the model. In 
section 3 we discuss the role of the Wilson mass parameter that enters the definition 
of the Dirac operator. In particular we show an interesting relation between the 
zero modes of the Dirac operator and the number of fermion flavors that the theory 
describes for background gauge field configurations. In section 4 we discuss the 
connection between topology and index of the Dirac operator. This was recently 
studied for semi-classical gauge field configurations in Ref. [E3]. Here we understand 
this connection as a function of the gauge coupling. We then argue that this is 
responsible for reproducing the singularities of the quenched approximation. In 
section 5 we calculate the spectral density of the Dirac operator as a function of the 
fermion mass. We also calculate the chiral condensate and the meson mass after 



5 A rather thorough analysis of the chiral condensate for the two flavor theory in the domain 
wall formalism has been done in Ref. fL7[ . 

6 Actions that have only finite range interactions are referred to as ultra-local. 



7 I thank Pavlos Vranas for reminding me of this fact. 
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including the fermion determinant and compare it with the results in the Schwinger 
model. A short discussion showing the relevance of the results to QCD and directions 
for future work are pointed out in section 6. 

2. The Model 

The model that we investigate can be described by the action 

5 = / 3 E( 1 -^[^ + 4])+E + m5 x J% (2) 

V V Z J x,y 

where Up is the phase obtained by multiplying the four £7(1) phase angles attached 
to a plaquette and (3 = 1/g 2 is the gauge coupling. Fermions of mass m interact 
with the gauge fields through the Dirac operator D, chosen to be the one proposed 
by Neuberger 0. In terms of the hermitian Wilson-Dirac operator 

H x , y = (6 XlV -nUl(x) [1 - 7m ] + 5 x , y+ My) [1 + 7 „]) + 5 x , yl5 {2 + M), (3) 

with the mass parameter M, D is given by 

D = l + jaH/VWH. (4) 

It is easy to show that D satisfies eq.([l|). The fermions obey (anti)periodic boundary 
conditions in the (temporal) spatial directions and the gauge fields are periodic. 

The condensate that characterizes chiral symmetry breaking was defined in Ref. 
|J. After integrating out the fermions this can be written as 

)/<[Det(£> + m)f') (5) 

where L s , L t are the lattice sizes in the spatial and temporal directions and Z takes 
care of the wave function normalization. Nf is the number of flavors. Nf = 
corresponds to the quenched theory. We do not consider Nf > 2 in this article. 

All our calculations are performed on an 8 x 16 lattice. Gauge fields are generated 
by a combination of head bath and over relaxation algorithm. Every gauge field 
configuration is generated after 500 sweeps (one sweep = one heat bath update + 5 
over- relaxation step over the entire lattice). The fermion determinants are absorbed 
into the observable and the expectation value is taken over gauge field configurations 
generated with the Boltzmann weight of the pure gauge theory. Although not a 
very reliable method for large lattices, the procedure seems to produce results with 
controlled errors in the present calculation. This is because the physics of the one 
flavor model is mainly sensitive to topology and the fermion determinant within 
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a topological sector does not fluctuate much. Further these sectors are updated 
well with the above algorithm on lattice sizes such as ours. Calculations involve a 
statistics over 2500 independent gauge field configurations. The errors are obtained 
through a jack-knife procedure. 

The model contains three additional parameters. The coupling g = 1/ yffi con- 
trols the lattice spacing. Since g is dimensionful in two dimensions, g ~ I /a. The 
continuum limit of this model is reached in the (3 —>■ oo limit. The mass parameter 
M does not effect the fermion mass, which is determined by m. It controls the 
effects of the fermion doubling. We will clarify its role further in the next section. 
Since it is difficult to control the four Fermi- couplings induced by the above model, 
it is considerably more difficult to do a systematic continuum extrapolation. Based 
on results in Ref. ]19j it is likely that in the continuum limit the above model will 
show differences with the results in the Schwinger model at the 10% level. Hence we 
study the model at a fixed (3, which also introduces finite lattice spacing errors. As 
we will see the results from the present study differ from the Schwinger model by 
10% — 15%, and based on the above discussion can be considered understandable. 



3. The role of M 



The parameter M is the mass of the Wilson Dirac operator that enters the 
definition of D. The theory can describe massless fermions for a range of values 
of M. Interestingly, as M changes it can change the number of massless fermions 
in the theory and thus can reintroduce fermion doubling. This happens at some 
critical values, say M c . It is easy to understand this in the free theory by analyzing 
the propagator. The free propagator at M — — 1 was constructed in Ref. [24] J 
Extending this result to arbitrary M we get 



_ 1 _ ^sin(^) 
Sf ^~2 2[N(p)+u(p)) (6) 



where — tt < p < tt, with 



<p) = 2 + M-£cos(p M ), N(p)= M^+^sin 2 ^) (7) 

When M > there are no massless particles. For —2 < M < the only pole in 
the propagator describing massless particles is at the origin of the (pi,P2) plane. 
However, in the region — 4 < M < — 2 two additional poles appear at (tt, 0) and 
(0, 7r). Here the Dirac operator describes three flavors of massless fermions. Finally 
for M < —4 a fourth pole at (n, tt) also appears and the theory describes four mass- 
less fermions. Thus the theory can describe different phases which are distinguished 

8 The free propagator for a whole class of Ginsparg- Wilson fermions was discussed in Ref. [B3[ . 
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by the number of massless fermion flavors. In the free theory the critical values 
M c = 0, —2 and —4 separate these phases. The number of massless flavors in the 
various phases discussed here differ from the those in domain wall fermions when the 



extra dimension is discrete [25]. This question was recently investigated by Shamir 



who observes the above results in the limit where the lattice spacing in the extra 



dimension becomes zero [26|. 



The critical values of M c will fluctuate in the presence of gauge fields. It is 
important to find a region in M in the full interacting theory where one can isolate 
the one flavor phase. In the weak coupling limit such a phase is most likely to 
exist. Since this also coincides with the continuum limit, both for the present model 
and in QCD, it is likely that one can actually obtain the physics of one flavor, 
sufficiently close to the continuum. On the other hand in realistic calculations it 
will be important to establish this before the results can be trusted. We will show 
that there is an interesting connection between the index (or equivalently the number 
of zero modes with a given chirality) and the number of massless fermions described 
by the Dirac operator in a background gauge field configuration at a given value of 
M. It may be possible to use this connection to determine the values of M that 
belong to the one flavor phase for each given background gauge field configuration. 

Consider a background gauge field configuration with a non-trivial index in the 
one flavor phase. As we cross a critical value of M the index can change. Such 



effects have been studied recently in Ref. [p7| . The change appears to be related 
to the change in the number of massless flavors the Dirac operator describes. We 
can predict it as a function of M based on this relationship. For example, let v ^ 
be the index of D in the one flavor phase. In the three flavor phase, we expect 
massless fermions from the poles at (0,0), (vr,0) and (0, 7r). However, the chirality 
definition of the fermions with poles at (tt, 0) and (0, tt) are opposite to that of (0, 0) 
and (tt, 7r). This means that in the three flavor phase the index of D will be —v. 
A similar argument shows that the index of the operator in the four flavor case 
would be zero. Since the operator is designed to produce zero modes equal to the 
index, this means that the Dirac operator will not have any (except accidental) zero 
modes in the zero and four flavor phases and will have v zero modes of opposite 
chirality in the one and three flavor phases. Thus one can scan through M to find 
the four phases for any gauge field with a non-zero index. In order to check the 
above observations numerically, we generated a series of gauge field configurations 
at (3 = 3. on a 8 x 16 lattice. In table [8] we give values of the index of the Dirac 
operator for ten such configurations for four values of M in the various phases. The 
results indicate that M — —1.0 is a region dominated by the physics of one fermion 
flavor, when /3 > 3. 

As far as we know, the above connection between the index of the Dirac opera- 
tor and the number of fermion flavors has not been appreciated before. In fact this 



6 





Index 




1V1 — l.U 


1V1 — -l.U 


1V1 — -o.u 


1V1 — -o.u 


U. 1UUO 


n 
U 


-i 


1 
1 


n 
U 


fl 1 

U. IOOj 


o 


1 


-1 


o 


0.1962 





3 


-3 





0.1669 














0.1822 














0.2309 





-1 


1 





0.1933 





2 


-2 





0.1904 














0.1753 





1 


-1 





0.1527 















Table 1: The table shows the index of the Dirac operator as a function of the Wilson 
mass M. The connection of the index with the number of massless fermion flavors 
and their chiralities is responsible for the above results. The first column indicates 
the value of the plaquette term of eq. (^) for the various configurations used. The 
results are consistent with theoretical predictions. 



connection leads to an amusing conjecture. In the region M < —4 there are four 
massless flavors but the Dirac operator does not produce any zero modes related to 
topology. The zero modes, if they exist, are accidental. However, there is a U(l) 
chiral symmetry in the action because the operator obeys the Ginsparg- Wilson rela- 
tion ||. Since there are no zero modes to break the symmetry through the measure, 
the exact chiral symmetry present must be flavored and hence not anomalous. This 
reminds one of staggered fermions where such a symmetry exists although the theory 
describes two flavors of massless fermions . 

In addition to producing the different phases, M is expected to renormalize the 
relevant couplings including the gauge coupling and the fermion mass. In fact it 
also introduces a wave-function renormalization. This can be seen by expanding the 
fermion propagator in small momenta near the poles in the various phases. We get 

!\M\{Pi + P2) + ••• expansion near (0, 0), M < 

\2+M\(pi + P\) + ••• expansion near (0, tt) or (tt, 0), M < -2 
u+m\ (Pi + P2) + ••• expansion near (tt, tt), M < -4 

(8) 

This shows that the wave-function renormalization at tree level turns out to be 
Z = 1/\M\ in the one flavor phase. As we will see below, the dominant effect 
of M on the chiral condensate of eq.@ comes from Z. Interestingly, we see that 
the residue at the various poles are different showing that the various massless 
fermions can have different wavefuntion renormalization constants. This suggests 
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that the renormalization properties of the three and four flavor phases needs to be 
understood further before they can be treated as a regularized version of legitimate 
field theories. 



4. Zero Modes and Quenching 



The one flavor Schwinger model is mainly governed by the dynamics of the gauge 
field topology and its connection with zero modes of the Dirac operator through the 
index theorem. Previous studies have shown that Ginsparg- Wilson fermions can 
reproduce the index theorem easily. We confirm this here. We choose the geometric 
definition of the topological charge given by 

Q = rE lo #p) (9) 

where the value of the log is defined on the principal sheet, i.e., < log(t/p) < 2tt. 
For smooth fields one expects that the index of the Dirac operator will be related 
to the above definition of the topological charge |24|. For rough fields there is no 



reason for such a relation. In order to quantify this effect we compared Q with 
the index of the Dirac operator at a strong coupling (3 = 0.5 and a moderately 
weak coupling /? = 3.0 lattice. The results for fifty configurations are plotted in 
the figure ([I]). These results show that even though the index of the Dirac operator 
is not related to the geometrical definition in any formal sense, the two start to 
agree quite well at moderately weak couplings and that (3 = 3 is sufficiently weak 
to reproduce the topological properties of the two dimensional theory. In figure (0) 
we plot the histogram for the topological sectors obtained from an ensemble of 2500 
configurations at (3 — 3. This further shows that the algorithm samples the various 
topological sectors adequately. 

The realization of the index theorem on the lattice has important consequences. 
In QCD it is well known that the mass of the rf particle is a consequence of the 
dynamics of the zero modes of the Dirac operator. In the quenched theory the 
zero modes are enhanced since the suppression of these modes due to the fermion 
determinant is absent. This causes many interesting singularities in the theory. In 
the infinite volume limit, these singularities have been predicted based on quenched 
chiral perturbation theory in Ref. [[$|. At small volumes semi-classical expansions 
suggest that there will be instanton like configurations that will dominate the path 
integral and this will lead to exact zero modes of the Dirac operator. This leads to 
new singularities. One of the consequences is the divergence of the quenched chiral 
condensate. 

The divergence of the quenched chiral condensate in the Schwinger model has 
been discussed in Ref. p9| . The condensate was first calculated in Ref. |3(| using 
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staggered fermions on a 16 x 16 lattice. The divergence as a function of the fermion 
mass was not found and was completely ignored even in the discussion. The singu- 
larity was first discussed in some detail in Ref. [[Jl[] where simulations on 64 x 64 
lattices were reported. These calculations have not been confirmed by other groups. 
Recent studies in QCD suggest that such singularities are difficult to obtain with 
staggered fermions f32| , [33|1 , but appear rather easily in the domain wall formalism 
13|| . The above discussion of the index theorem shows that with Ginsparg- Wilson 
fermions these divergences must appear easily. In order to show the dramatic dif- 
ference between staggered fermions and the Ginsparg- Wilson fermions, we plot the 
quenched chiral condensate as a function of the fermion mass in figure (|3|) for both 
types fermions. 



5. Dirac Spectral Density and Meson Mass 



The study of the chiral condensate and the meson mass is important in any study 
of the Schwinger model. The chiral condensate is given by eq.([5|). Written in terms 
of the eigenvalue density p(9) of the Dirac operator |J, it reads 

, tm ^ r ^ / n\ 1 + cos(0) + ism(9) . . 

(ibib) = Z d9 o(e)— — ^ — ^ r. 10 

7-Tr V ; 2(1 -cos(0) -isin(0) + am) v ' 

In deriving the above representation, we have used the fact that the eigenvalues of 
the Ginsparg- Wilson Dirac operator^ are distributed on a unit circle in the complex 
plane centered at (1,0). The density is normalized using the relation J d9p(9) = 2. 

An eigenvalue at 9 = is an exact zero mode. It is easy to see that every 
eigenvalue at 9 = will be accompanied by an eigenvalue at 9 = tt for topological 
reasons. Separating these topological zero modes, we find 

p(9) = p [6(9) + 6(9 - tt)] + p\9) (11) 

In figure (f|) we plot p'(9) at (3 — 3.0 and \M\ = 1.0 obtained using quenched 
configurations. For comparison we also plot the free field eigenvalue density. 

It is useful to understand the dependence of p(9) on the dynamical fermion mass. 
Since all the eigenvalues are measured, this is easily obtained. For each configuration 
we divide the range (— n, n) into 50 bins and count the number of eigenvalues falling 
into each bin and divide it by (L s L t )2n/50 to obtain the value of p'(9) where 9 
refers to the center of the bin. For 9 = we give the value of po which is obtained 
by counting the zero modes and dividing by the L s L t . p(9) for four values of the 
dynamical fermion mass is given in table |. The errors are obtained by a jackknife 
method. p(9) is even and we give the values for positive 9. From the results it 

9 We assume that D obeys = 75D75. 
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appears that the statistics of 2500 configurations is far from sufficient to see the 
effects of the fermion determinant except in the region \6\ < 0.35. This corresponds 
to masses less than 0.06 in lattice units. 

The chiral condensate is evaluated for three different values of M, in the one 
flavor phase, in order to study its effects. The dominant effect comes from the 
tree level wave function renormalization factor Z = 1/\M\. This factor will get 
corrections from loop effects. The results after including this factor are shown as a 
function of the fermion mass in figure (^). In the continuum it can be shown that 

where C = 1.78107... is the Euler's constant. This is indicated on the m = axis. 
As discussed in section 2, it will be difficult to reproduce the results of the Schwinger 
model with errors of less than 10 — 15% from the present lattice study. Our results 
are consistent with this expectation. 

Next we look at the two point meson correlator as defined by 
G„(t) = ^E((tr{(^0| 75 ^|,,t)(,,t| 75 ^l-|,,0)} 

- tr{<x ) 0| 7 b^^|a: ) 0>}tr{(y ) t|7b^ ; |y ) «>}) [Det(£> + m)]"') 
/([Det{D + m)) N f) (13) 



in order to extract the meson mass fi. The kets \x,t) represent the normalized 
position basis states and the trace is over the spinor index. Using the eigenvalues and 
eigenstates of D it is easy to evaluate the propagators necessary for the calculation 
of G nn (i). In figure (§) we plot it as a function of t. The function appears to fit to 
the usual form of v4cosh(m(t — T/2)), suggesting that conventional methods will be 
useful to extract the meson masses. A simple estimate of the meson mass can be 
made using the effective mass analysis. Defining the effective mass by 

/iefr(i) = - \og(G wn (t)/G nw (t - 1)). (14) 

we plot the results in figure (|7|), which shows evidence for plateaux. We chose the 
values of fi e s at t — 5 to estimate the meson mass. The meson masses are plotted in 
figure (|S|) as a function of the fermion mass. The continuum result in the Schwinger 
model, fi/g = l/y/ft is shown on the m = axis. Again our results agree with in the 
expected 10 — 15% error margin. It is interesting to note that no factors of Z are 
necessary here as this is a spectral quantity. However, the dependence of the meson 
mass on M indicates a renormalization of the fermion mass and the coupling. This 
effect can be studied more concretely and is left for the future. 
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6. Discussion and Conclusions 



We have reproduced the essential features of the Schwinger model rather easily 
on a coarse lattice using Neuberger's formulation of Ginsparg- Wilson fermions. This 
was possible because the physics of the one flavor model is essentially governed by 
the topological sectors, which was reproduced quite well on a moderately coarse and 
finite lattice. We discovered an interesting connection between the zero modes and 
number of fermion flavors as a function of M. The value of M was important in 
getting the right connection between the number of flavors and the zero modes of 
the Dirac operator. In the region > M > —2, where the theory is expected to 
describe the right physics the dominant effect of M in the chiral condensate was 
a wave-function renormalization. Evidence for the renormalization of the fermion 
mass and the coupling was also seen. 

All this means that questions that are dominated by topological effects in QCD, 
such as the rj mass, may also be easily addressed with these fermions, once the right 
region in M is isolated. Reproducing the physics of the anomaly has other interesting 
consequences. In particular it can help in a better understanding of the QCD phase 
transition. Most calculations in lattice QCD indicate a second order phase transition 
for two flavor QCD. The most reliable calculations come from staggered fermions 
which have a remnant chiral symmetry. However, we have seen that staggered 
fermions produce very weak anomalous effects. These results have been confirmed in 



QCD just above the phase transition [|34| . As is well known |35j , a very weak anomaly 
close to the transition can make the transition first order. Since this is not seen in 
simulations, our understanding of the two flavor transition may not be complete. 
On the other hand, since the Ginsparg- Wilson fermions reproduce the effects of 
the anomaly rather easily and in addition have the correct flavor structure, they 
may give a more definitive understanding of the chiral transition. First such results 
have been reported in the domain wall formulation j36]. The present study also 



indicates that the Ginsparg- Wilson fermions behave quite differently with respect 
to quenching: another feature of reproducing the zero mode physics well. Hence 
the effects of the fermion determinant will play a more interesting role in QCD and 



quenched singularities |23 could be easily seen with these fermions. This is perhaps 



the easiest to study, since it involves only quenched calculations. 

One of the striking feature of QCD is the spontaneous breaking of chiral sym- 
metry. Unfortunately this cannot be addressed in a two dimensional model due to 
the Mermin- Wagner theorem. On the other hand, there is very interesting physics 
associated with the many flavor Schwinger model. For example one expects a non- 
analytic dependence of the meson mass and the chiral condensate on the fermion 
mass |2^| . As far as we know this has not been reproduced in the past mainly due to 
the lack of a good fermion formulation^ We plan to extend the present calculation 

10 This has been discussed for the first time with domain wall fermions in Ref. | fl7| 
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in that direction. However, since we cannot expect non-analytic behavior to emerge 
on a finite lattice, a more careful finite size scaling analysis would be necessary. Thus 
the many flavor Schwinger model will be considerably more difficult, but perhaps 
can serve as a useful guide for the study of chiral symmetry in QCD. 
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Table 2: The table gives the eigenvalue density of the Dirac operator as a function 
of the fermion mass that appears in the determinant. The last column refers to the 
infinite mass limit (which is equivalent to the quenched limit) and is shown in figure 
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Figure 1: The geometric definition of the topological charge is compared with the 
index of the Dirac operator for 50 configurations generated at ft = 0.5 (top) (3 = 3.0 
(bottom). 
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Figure 2: This figure shows a histogram of the topological sectors generated at (3 = 3.0 
in the 2500 configurations that were analyzed. 
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Figure 3: TTie chiral condensate including the factor Z = 1/\M\ in the quenched 
limit is shown. For comparison we also include the results from staggered fermion 
simulations. 
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Figure 4: This figure shows the eigenvalue density of D for the quenched run at 
(3 = 3.0 and for free field theory. 
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Figure 5: The chiral condensate including the factor Z = 1/\M\ is plotted as a 
function of the fermion mass for various values of M. The dominant role of \M\ 
comes through Z . The continuum value in the massless Schwinger model is shown 
on the m/g = axis. 
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Figure 6: 77ms figure shows the zero momentum correlation function G n7T (t) as a 
function of time slice. 
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Figure 7: The effective meson mass at various time slices is plotted here to show 
evidence for plateaus. 
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Figure 8: The meson mass is plotted as a function of the fermion mass for various 
values of M. The value in the continuum massless Schwinger model is represented 
by a star at m/g = 0. 
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